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1.  Summary 


->(>>> 


This  report  summarizes  work  done-  by  the  authors  during  the  past  calendar  year  on 
the  dynamics  and  control  of  flexible  and  articulated  spacecraft.  The  combined  dynamical 


effects  of  elasticity  and  a  rotating  reference  frame  have  been  explored  for  structure  in  a  zero 


gravity  environment.  A  simple  yet  genera]  approach  to  modeling  has  bee^ developed,  and 

Jhis  approaeh-has-beefpapplied  to  analyze  <uv  detaifnhe  dynamics  of  a  specific  prototypical 

Oju  *  ;  r 

structure.  The  effects  of  energy  dissipation  fitave -been  included  and  studied  in  depth  for  a 

v...  l"  I  ' 

model  problem.  Equilibria,  bifurcations,  and  asymptotic  stabilitjyhave-beei^analyzed  in 

some  carefully  chosen  examples  which  capture  the  essential  general  features  of  nonlinear 

I  / 

distributed  parameter  models  of  rotating  elastic  structures.  '  X  ■ .  '  J 


2.  Research  Objectives 


The  main  focus  of  this  research  program  has  been  on 

•  Development  of  a  qualitative  theory  of  the  dynamics  of  rotating  elastic  and  articu¬ 
lated  structures; 

•  Development  of  simplified  models  and  an  understanding  of  their  domains  of  validity; 

•  Development  of  functional  analytic  techniques  for  analyzing  the  evolution  of  the 
relevant  distributed  parameter  systems;  and 

•  Development  of  an  appropriate  control  theory  for  nonlinear  distributed  parameter 
systems  of  the  type  required  to  model  rotating  elastic  structures. 

A  general  approach  to  modeling  rotating  elastic  structures  has  been  developed  and  is 
described  in  subsection  3.2.  Here  it  is  shown  that  mechanical  effects  due  to  the  rotation  of 
our  reference  frame  may  be  modeled  independently  of  affects  due  to  deformations  of  the 
elastic  material.  In  subsection  3.4  we  present  what  we  believe  is  the  simplest  reasonable 
continuum  mechanical  model  of  a  beam  undergoing  three  degree  of  freedom  rotations. 
Subsection  3.5  and  3.6  treat  stability  issues  and  teh  steady  state  behavior  of  such  models. 

Future  research  will  focus  on  development  of  techniques  for  actively  controlling  the 
motions  of  rotating  structures.  In  particular  it  is  our  aim  to  understand  how  commanded 
changes  in  angular  momentum  can  be  used  to  steer  a  rotating  elastic  structure  from  one 
equilibrium  state  to  another.  The  relevant  control  theory  of  nonlinear  distributed  param¬ 
eter  systems  will  be  developed  against  the  backdrop  of  this  specific  problem. 


3.  Status  of  Research 


3.1.  Introduction 

There  "is  now  a  fairly  general  awareness  among  aerospace  engineers  that  the  dynamics 
of  the  coifiplex  spacecraft  currently  in  production  and  on  the  drawing  boards  will  be 
greatly  influenced  by  continuum  mechanical  effects  such  as  elasticity.  Indeed  as  the  designs 
being  contemplated  increase  in  size  and  complexity  ([13])  the  dynamic  effects  of  flexible 
members  become  more  important,  and  recognizing  this,  many  researchers  over  the  last 
decade  have  focused  their  efforts  on  obtaining  new  methods  for  the  design,  analysis,  and 
control  of  flexible  mechanisms.  Space  does  not  permit  (nor  would  it  be  in  keeping  with 
the  main  purpose  of  this  report)  to  survey  the  vast  literature  on  flexible  space  structures; 
the  interested  reader  can  get  some  idea  of  research  activity  in  this  area  by  referring  to  any 
one  of  a  number  of  collections  of  papers  and  conference  proceedings,  such  as  [12]. 

The  purpose  of  this  section  is  to  describe  recent  research  which  has  been  aimed  at 
developing  a  mathematical  theory  of  the  rotational  dynamics  of  complex  mechanical  sys¬ 
tems  which  include  articulated  and  elastic  components.  Our  objective  in  this  research  has 
been  to  carry  out  a  study  of  the  global  qualitative  dynamics  of  such  systems  in  sufficient 
depth  as  to  allow  predictions  regarding  the  stability  and  asymptotic  behavior  of  spacecraft 
due  to  a  variety  of  energj  dissipation  mechanisms  such  as  viscoelastic  material  damping 
of  vibrations  of  elastic  parts.  We  believe  that  historical  evidence  points  to  the  value  of 
developing  a  fairly  complete  global  asymptotic  stability  theory  of  this  type  since  there 
are  numerous  examples  of  missions  in  space  which  did  not  achieve  their  stated  objectives 
because  certain  long  term  mechanical  effects  were  never  adequately  taken  into  account  in 
the  mission  planning.  Explorer  I,  the  first  successfully  launched  American  satellite,  pro¬ 
vides  the  best  known  example  of  such  untoward  behavior.  Upon  achieving  earth  orbit, 
the  pencil-shaped  satellite  was  supposed  to  rotate  about  its  major  axis  of  inertia.  Before 
it  had  completed  one  orbital  revolution,  however,  radio  signals  indicated  that  a  tumbling 
motion  had  developed  and  was  increasing  in  amplitude. 

While  explanations  of  Explorer’s  errant  behavior  have  been  offered  by  a  number  of 
researchers  (see  e.g.  references  to  this  problem  in  [7]),  we  are  aware  of  no  attemps  to  obtain 


a  rigorous  mathematical  analysis  of  such  occurences,  with  the  exception  of  [3],  where  some 
of  the  results  of  the  present  paper  were  announced.  We  also  mention  a  recent  paper  [8] 
which  deals  with  a  somewhat  different,  although  related  aspect  of  a  similar  problem. 

The  present  section  is  organized  as  follows.  In  subsection  3.2  we  derive  equations 
describing  The  rotational  dynamics  of  complex  structures.  These  include  the  general  effects 
of  inertial  forces  created  by  rotation  of  the  reference  frame.  Subsection  3.3  focuses  on  a 
general  theoretical  framework  for  Lagrangian  mechanics  with  damping.  In  subsection  3.4, 
a  simple  structure  consisting  of  a  rigid  body  with  an  elastic  beam  appendage  is  studied, 
and  we  present  what  we  believe  is  the  simplest  reasonable  continuum  mechanical  model  of 
such  a  system  undergoing  three  degree  of  freedom  rotations.  The  asymptotic  steady  state 
dynamics  for  this  system  are  studied  in  subsection  3.5,  and  in  subsection  3.6  we  present 
a  detailed  analysis  of  the  qualitative  dynamics  of  a  closely  related  model  having  only  one 
rotational  degree  of  freedom. 


I  lit  1 


8.2.  The  rotational  dynamics  of  complex  structures. 

In  this  section  we  will  derive  equations  of  motion  for  a  class  of  structures  consisting 
of  elastic,  fluid  or  rigid  components.  While  these  equations  are  completely  general,  they 
are  most  useful  in  describing  any  structure  whose  configuration  is  conveniently  specified 
by  the  position  of  the  structure  with  respect  to  a  moving  coordinate  frame  together  with 
the  position  of  that  frame  in  space,  i.e.  relative  to  some  fixed  inertial  system. 


Figure  2.1:  The  body  frame  is  translated  und 


rotated  with  respect  to  the  inertial  frame. 

We  fix  an  orthonormal  basis  forming  an  inertial  frame  (“space  frame”),  and  choose  a 
“body”  coordinate  system  designated  by  a  set  of  orthonormal  vectors  gi.g2,g3-  Choice 
of  the  g,’s  will  depend  on  the  particular  problem  at  hand,  and  it  affects  the  simplicity  of 
resulting  equations.  In  subsection  3.4,  where  we  treat  an  example  of  a  rigid  body  with  an 
elastic  appendage,  we  affix  the  body  frame  g,  to  the  rigid  body,  although  other  choices  are 
possible  [6] . 

The  position  and  orientation  of  the  body  frame  (with  respect  to  the  chosen  inertial 
frame)  may  be  described  at  each  time  t  by  an  element  of  the  special  Euclidean  group, 
5£'(3,R),  represented  by  a  4x4  matrix 

x(0-(T  «f) 

(see  figure  2.1),  where  V(f)  £  50(3)  is  an  orthogonal  matrix  describing  the  orientation  of 
the  body  frame  and  y  £  R3  is  the  position  of  the  origin  of  that  frame  in  space,  (cf.  [5]). 

Therefore,  if  a  point  P  is  given  by  the  vector  a  in  the  body  frame  and  by  vector  U  in 


the  space  frame,  then  u  and  U  are  related  via 


U  =  Xu, 

where  U  =  (Y)  and  u  —  (“). 

The  positions  of  various  elements  of  our  flexible  structure  relative  to  the  body  frame  T  - 
will  be  described  by  a  vector  function  u(z,t)  =  u(z,,z2,z3,t)  denoting  the  position  at  time 
t  of  each  particle  whose  “unperturbed”  position  is  at  z.  Here  “unperturbed”  can  mean 
either  initial  or  undeformed,  depending  on  one’s  choice.  In  subsection  3.4,  z  will  denote 
the  neutral  position  of  a  particle  of  the  flexible  structure. 

In  summary,  we  consider  systems  whose  configuration  space  is  given  by 
{<?}  =  {(1.  u,  y)}  =  SE(3)  x  C  =  50(3)  x  R3  x  C,  where  C  =  { !/(•,<)}  is  a  suitably  defined 
function  space  whose  elements  are  functions  u(z)  describing  the  configuration  of  the  body 
relative  to  the  body  frame. 

We  describe  now  the  kinematics  of  the  system  and  give  an  expression  for  the  kinetic 
energy.  The  evolution  of  the  matrix  A"(t)  €  SE[ 3)  can  be  described  by  a  differential 
equation 


X(t)  = 


V(<)n(<)  y(<) 


0 


where  fl(/)  is  the  skew-symmetric  matrix 


(0  -w3  u>2  \ 

W3  0  —  Wj  I 

—  u>2  U/j  0  J 

of  angular  velocities  about  the  corresponding  body  axes  g,,  i  =  1,2,3. 

The  inertial  coordinates  U(z,t)  of  a  point  P  are  related  to  its  body  coordinates  u(z,t) 
via 

U(z,l)  =  Xu{z,l) 

and  the  corresponding  velocities  are  given  by 


jU  =  Xu  +  Xu,  = 

_  (  Y(Clu  +  ut)  +  y 


=  ( 


)■ 


where  ut  denotes  the  partial  derivative  with  respect  to  t. 
The  kinetic  energy  of  the  body  is  then  given  by 


T(l,g)  =  \JjUH2dm=lJBt\Y(nu  +  u,)  +  yl\Um,  (2.1) 

where  B  denotes  the  point  set  comprising  the  body  (at  time  t)  described  in  the  body 
coordinate  system,  and  dm  is  the  mass  distribution  in  this  coordinate  system. 

The  kinetic  energy  of  almost  any  rotating  structure  will  be  of  this  form,  as  it  does  not 
depend  on  the  constitutive  relations  governing  the  structure  itself.  More  refined  dynam¬ 
ical  models,  as  treated  in  subsequent  sections,  will  embody  structural  information  in  the 
expression  for  potential  energy 

V(q)=V(Y,u,y) 

(which  we  assume  to  be  independent  of  g).  Without  specifying  the  form  of  V  at  this  point, 
we  prove  the  main  theorem  of  this  section. 

Theorem  2.1.  Let  A  be  defined  by 


A{t)  =  I  ((Hu  +  ut)uT)adm, 
J  B 


where  Ma  =  \{M  —  MT)  denotes  the  anti-symmetric  part  of  a  matrix  M.  Equations  of 
motion  of  any  system  whose  configuration  space  is  {q}  =  {(T,  u,y)}  =  50(3)  X  C  X  R3  as 


above  are  given  by 


Ml)  +  [n(<M(01  -  l(YTyuT)adm  =  T  -  (Y-'V Y)c 
J  B 

utt  +  n2u  -f  flu  +  2Hu(  4-  YTy  =  F  — 

/dV 
Y{Q2u  +  flu  +  2Hu|  +  utt)  +  ydm  =  f  - 


y  =  Yn 


(2.4)' 


where  Vy  is  the  matrix  of  partial  derivatives  ,  and  jL  {s  the  Frechet  derivative  of  V 
with  respect  to  u. 


In  the  skew-symmetric  matrix 


T  = 


/  0  -r3  t2  \ 

t3  0  -r, 

V  -T2  r}  0,  J 


t,  is  the  net  nonconservative  torque  applied  to  the  system  about  the  body  axis  g,,  F  =  .F(2) 
is  the  distributed  nonconservative  force  density  acting  on  the  particle  positioned  at  u(z,t) 
expressed  in  the  body  coordinate  system,  and  f  is  the  net  nonconservative  exogenous  force. 

Remark:  In  subsequent  sections  nonconservative  forces  will  arise  due  to  viscoelastic 
damping. 


Before  giving  the  proof,  we  rewrite  these  equations  so  as  to  provide  a  clearer  picture 
of  the  physical  situation.  Let  w  =  (uq  ,u>2,uj3)t  be  the  angular  velocity  of  the  body  frame 
expressed  in  that  frame-,  u;  is  related  to  the  angular  velocity  matrix  fl  =  Y~'Y  as  follows: 


l  0 

—  W3 

u.>2  ^ 

^3 

0 

-UJX 

,  w  = 

u2 

V  “^2 

W1 

0  ) 

{^3/ 

Denote  by  5  the  operator  taking  a  skew-symmetric  matrix  D  into  vector  u>  :  Sf)  =  jj.  One 
easily  verifies  the  following: 

Lemma  2.1.  Given  any  pair  A,B  of  skew- symmetric  3x3  matrices  and  any  pair  u,v 
of  3-vectors,  the  following  identities  hold 


(i)  S(\A.B])  =  S(A)xS(B) 

(ii)  S(uvT  —  vuT)  —  1  xu 

(iii)  Au  =  S(A)  x  u 

Applying  the  operator  5  to  equation  (2.2)  and  using  Lemma  2.1  on  equations  (2.2)-(2.4) 
we  obtain 


Corollary  2.1  Equations  of  motion  (2.2)-(2.4)  are  equivalent  to 

a{t)  +  u(t)  x  a[t)  +  f  u(z,t)  x  Y~1y(t)dm  =  5(1  -  (>' _  1  V”v )a )  (2-5) 

J  B 

9 


sv 


Utt  +  u  x  (*'  x  u)  +  ij  x  u  +  2ui  x  Uj  +  y  1  y  =  F  -  j- 

f  dV 

J^[Y(u  x  (u  x  u)  +  (j  x  u  +  2^'  x  u(  +  utt)  +  y\dm  =  /  - 

u’here  a  is  'defined  by 

.**■  a(t)  =  /  t<  x  (tq  +  x  u)dm. 

J  B 


(2.6) 

(2.7) 


Written  in  this  way  equations  (2.5)-(2.7)  give  an  explicit  description  of  the  inertial 
forces  on  the  mechanical  system  viewed  in  the  moving  body  frame.  Introducing  the  deriva¬ 
tion  D  —  4  (•)  -I-  u  x  (•),  we  obtain  yet  another  rendering  of  these  equations. 


Corollary  2.2  Equation  of  motion  ( 2.2)-(S .4)  o.re  equivalent  to 

f  ux  (D2u  +  Y~'y)dm  =  S{7  -  {Y-'Vy)a) 
J  B 


D*u  +  Y-'y=  F  -  ^ 


6V 


6u 


r  QY 

JB(YD2u  +  y)dm  =  f  - 


(2.5) ' 

(2.6) ' 

(2.7)' 


The  motions  of  any  complex  structure  undergoing  free  or  forced  rotation  are  described 
by  equations  (2.5)-(2.7).  This  formulation  is  thus  fairly  general,  and  it  can  incorporate 
external  forces  and  torques  (due,  for  example,  to  gravitational  and  magnetic  fields)  and 
internal  forces  (due,  say,  to  actuation  of  joints  or  the  constitutive  properties  of  the  mate¬ 
rial).  In  the  next  subsection  3.3  we  will  incorporate  dissipative  effects  in  this  formulation, 
and  in  subsection  3.4.  we  shall  develop  a  complete  dynamical  model  of  a  system  where  the 
constitutive  relations  are  those  of  a  simple  damped  beam. 


Proof  of  theorem  2.1:  The  proof  will  be  given  in  two  parts:  First,  the  treatment 
of  the  inclusion  V'  €  50(3)  as  the  holonomic  constraint  onto  50(3)  as  a  submanifold 
of  GL{ 3),  with  an  appropriate  modification  of  the  Lagrange  equations  (Lemma  2.2).  and 
second,  the  utilization  of  the  left-invariance  of  (the  rigid  path  of)  kinetic  energy  to  simplify 
the  resulting  equations  (Lemma  2.4). 

We  would  like  to  write  the  equations  of  motion  of  the  structure  in  the  Lagrangian  form 

d  dT  dT  dV  ...  . 

dt  dq  dq  ~  dq'  q~ 


10 


An  appropriate  modification  of  the  Lagrange  equations  is  needed,  however,  to  account  for 
the  fact  that  Y  is  constrained  to  the  submanifold  50(3)  C  G7(3).  Such  a  modification 
is  unnecessary  if  the  Lagrangian  equations  are  expressed  in  terms  of  a  local  coordinate 
system  on  the  constraint  manifold-cf.  [1],  page  77.  This  approach  ignores  the  symmetry 
in  our  system  of  equations,  however.  The  following  Lemmas  2.2  and  2.4  are  key  to  our 
proof  of  Theorem  2.1 


Lemma  2.2  Any  extremal  q(t)  of  the  action  f  L(q,q)dt  with  q  constrained  to  a  sub¬ 
manifold  A/q  of  a  Riemannian  manifold  M  satisfies  the  differential  equation 


d^dL 
M*  dq 


dL 


dq 


)=°, 


where  is  the  normal  projection  from  TqM  onto  T^Mq. 

We  omit  a  straightforward  proof.  Corresponding  to  the  three  components  of  q  = 
(Y.u.y)  we  obtain  three  components  for  the  equations  of  motion: 


PY  7F7 


d 

dT 

dT 

dV\  „ 

(2.8) 

dt 

dY  ' 

~  Fy  + 

dYj  ~  °’ 
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6V 

(2.9) 

dt 
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Su 

Su  ’ 
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dT 

dT  _ 

dV 

(2.10) 

dt 

dy 

dy  ~ 

dy' 

where  Py  is  the  orthogonal  projection  from  TyGl(3)  onto  TySO(3)  in  the  trace  norm 
( A,B )  =  trATB.  Here  .  |y  are  the  derivatives  with  respect  to  the  standard  Riemannian 
structure  on  TGL( 3):  they  can  be  represented  as  matrices  of  partial  derivatives:  ^  = 
and  |y  =  (^~).  denotes  the  Frechet  derivative  of  T  with  respect  to  the 
distributed  parameter  u.  The  expression  for  the  projection  Py  is  provided  by 


Lemma  2.3.  For  any  Y  £  50(3)  and  A  C  TyGl{3) 

PyA  =  Y(Y~'A)a, 

where  Xa  =  ^(X  -  XT) 


(2.11) 


is  the  anti-symmetric  part  of  X. 


Proof.  This  follows  from  decomposing  the  Lie  algebra  of  gl( 3)  as  the  orthogonal  direct 
sum  of  symmetric  and  skew  symmetric  matrices.  (gl( 3)  =  so( 3)  ©  soft)1 .)  | 


Leinnja  2.4  Suppose  T  is  a  left-invariant  function  on  TGl(n),  i.e.  suppose  there  is  a 
function  Kf^defined  on  gl(n)  (=  the  Lie  algebra  of  Gl(n)  =  space  of  real  n  x  n  matrices) 
such  that  r(Y,F)  =  K(Cl)  where  fl  =  Y~'Y.  Then  for  (Y,Y)  €  TSO{n) 


where  M  —  (^)a>  and  is  the  derivative  of  K  with  respect  to  fi  evaluated  at  fi  =  Y  IY. 

Remark.  There  is  an  orthogonal  direct  sum  decomposition:  gl[n)  =  so(n)  ©  so(n)-1-, 
where  so(n)  =  the  Lie  algebra  of  n  x  n  skew-symmetric  metrices,  and  where  orthogonality 
is  defined  in  terms  of  the  trace  form  inner  product  (A,  B)  =  tr  ABT  defined  on  gl(n).  If  the 
function  K  appearing  in  the  statement  of  Lemma  2.4  can  be  decomposed  as  K  =  K\  -fi 
where  K,  depends  only  on  the  i-th  component  (i  =  1,2)  in  this  orthogonal  direct  sum, 
then  M  =  K[{ fi),  where  by  K[  we  mean  the  derivative  of  A'j  with  respect  to  the  natural 
differentiable  structure  on  so(n)  defined  in  terms  of  the  Killing  form. 


Proof.  and  may  be  thought  of  as  elements  in  the  cotangent  bundle  T*Gl{n). 

Making  the  usual  identifications,  the  standard  Riemannian  structure  on  Gl(n)  may  be 
prescribed  explictly  in  terms  of  the  trace  form,  and  we  may  write  an<^  In  45 

n  x  n  matrices 


£dT 
dt  qy 
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For  (y,  Y)  €  T  SO(n),  we  have  (V  1)T  =  Y  and  nT  =  -fi,  and  the  result  follows  from 
Lemma  2.3.  | 


The  proof  of  Theorem  2.1  proceeds  as  follows.  Using  left-invariance  of  the  first  term  in 
the  expression  for  kinetic  energy  (cf.  (2.1)) 

T  =  ^  ||Hu  +  it/ )|2  -fi  2(fiu  -fi  U/ .  y - 1  y)  ||y)|2  dm, 

we  obtain  equation  (2.2)  as  a  consequence  of  Lemma  2.4.  The  remaining  equations  (2.3) 
and  (2.4)  follow  by  direct  computation  from  equation  (2.9)  and  (2.10).  We  omit  the  details. 


3.3.  Lagrangian  Mechanics  with  Damping 


A  major  advantage  of  Lagrangian  versus  Newtonian  mechanics  is  the  invariance  of 
Lagrange’s  equations  with  respect  to  coordinate  changes;  it  is  this  invariance  that  facilitates 
significantly  the  derivation  of  our  equations  of  a  motion.  It  is  thus  desirable  to  have  the  ’• 
extension  to  the  dissipative  case.  Such  a  modification  is  described  in  [9],  [l l] ;  we  reproduce  ' 
it  here  in  a  slightly  more  general  form. 

Let  D(q,q)  be  the  so-called  dissipation  function  defined  as  follows: 

qDq  =  rate  of  dissipation  of  energy  per  second. 

(One  can  think  of  Dq  as  the  generalized  dissipation  force,  and  q  =velocity.  The  above 
simply  says:  velocity  ■  force  =  power). 

Let  L{q,q)  be  the  Lagrangian  of  the  system.  The  equations  governing  the  system  are 

^  Lq  —  Lq  -f  Dq  =  0.  (3.1) 

Remark  3.1.  If  D  is  quadratic  in  q  (as  will  be  the  case  in  the  application  presented 
in  the  next  section),  then 

D  =  ^qDq  =  ^(rate  of  dissipation). 

Equations  (3.1)  are  consistent  with  the  definition  of  D ,  as  shown  by 
Theorem  3.1.  If  E(q,q)  is  the  total  energy  of  the  system,  then 

=  -qDq. 

Proof.  E  is  given  by  the  Legendre  transform*  in  q  of  L: 

r  dL  t 
E  =  ^-L- 

Differentiation  by  time  gives 

■  .dL  .d  ,dL .  _  .  _  . 

E  =  "si +  >  -  L'q  -  L*q  = 

*  Usually  E  is  expressed  as  a  function  of  q,  p  =  we  keep  q  rather  than  p  here. 


we  used  equation  (3.1)  in  the  last  step.  | 

Equation  (3.1)  has  the  same  invariance  property  as  the  conservative  Lagrange  equa¬ 
tions. 

Theorem  3.2. [  1 1]  The  dissipative  Lagrangian  system  (8.1)  is  invariant  under  the 
change  of  variables  q  =  q(Q).  More  precisely,  if  q(t)  satisfies  (8.1),  then  Q(t)  satisfies 
equations  of  the  same  form: 

jtSL( j-L<?  +  D(j=0)  (3.1)' 

where  L{Q,Q)  =  L(q(Q),  q'(Q)Q)  and  V(Q,Q )  =  D(q(Q),q'  (Q)Q). 

Proof:  A  simple  calculation  shows: 

■ 


3.4.  A  Rotating  Rigid  Body  with  a  Beam  Attachment 

Consider  the  spacecraft  depicted  in  Figure  4.1.  The  key  features  of  this  structure  are 
a  rigid  body  to  which  a  flexible  cantilevered  beam- like  appendage  of  length  £  is  attached. 

* 


Figure  4.1:  A  rigid  body  with  cantilevered  beam  attachment. 

As  we  have  mentioned  in  subsection  3.2,  we  affix  the  moving  frame  to  the  rigid  body. 
More  specifically,  we  place  the  origin  of  the  frame  at  the  point  of  attachment  of  the  beam 
and  align  the  z3  =  2-axis  along  the  undeflected  beam.  Viewing  this  cantilevered  beam  as 
essentially  a  one-dimensional  object,  we  describe  the  elastic  deformations  u(z3,t)  =  u(2.<) 
with  respect  to  the  coordinate  axes  (21,22,23)  =  (x,y,z)  depicted  in  Figure  4.1.  More 
precisely,  u[z,t)  is  the  position  of  the  particle  whose  neutral  position  is  at  (0,0,2).  The 
decomposition  of  the  system  into  the  rigid  part  and  the  elastic  beam  corresponds  to  the 
decomposition  of  kinetic  energy  (2.1)  into  the  sum  of  rotational  and  translational  energies 
of  the  rigid  part  and  the  energy  of  the  beam  (Note  that  u(-,f)  restricted  to  the  rigid 
component  is  just  the  identity  mapping  for  all  t .)  We  have 


where  S(fi)  —  w,  c  is  the  center  of  mass  of  the  rigid  body  in  the  body  frame,  mj,  is  the 
mass  of  the  rigid  body,  we  have  scaled  the  linear  mass  density  of  the  beam  to  be  one,  and 
the  inertia  tensor  with  respect  to  the  body  frame  is  given  by 

(h  hy  Izz' 

J  =  Izy  ly  Iyz  , 

i.  ^zz  lyz  Iz  , 

where  Iz  is  the  moment  of  inertia  with  respect  to  the  i-axis,  etc.  (cf.  [l]). 

Since  the  beam  is  clamped  at  the  origin  of  the  (*,  y,  z)-coordinate  system  and  free  at 
its  other  end,  the  following  boundary  conditions  are  assumed: 

U«(M  =  ^(<M)  =  (M)  =  ^r(U)  =  o,  i  =  1,2, 

and  u3(0,<)  =  ^(£,  t)  =  0-  These  boundary  conditions  are  standard  in  the  theory  of 
clamped-free  beams.  Here  ui,u2,u3  are  the  deflections:  u  =  (ui,u2,2  +  tt3).  Note  that  u3 
is  not  the  2-coordinate  of  u. 

The  equations  of  motion  for  our  rotating  satellite  are  obtained  by  incorporating  (2.5)- 
(2.7)  into  the  formalism  of  Lagrangian  mechanics,  as  discussed  in  the  previous  section. 
Thus  we  look  for  extremals  of  the  Lagrangian  L  =  T  —  V  with  kinetic  energy  T  given 
above  and  potential  V  given  by  the  strain  energy 


V(«)  =  \  /0'|«1  («?)’  +  '=§-z,  (4.2) 

where  only  quadratic  terms  were  retained  and  the  material  is  assumed  to  obey  a  linear 
Hooke’s  law.  Here  /xj  (respectively  /x2)  gives  the  bending  elasticity  within  the  12-plane  (yz- 
plane  respectively),  and  /x3  is  the  Hooke’s  constant  giving  the  beam's  stretching  elasticity. 
Unless  the  beam  is  abnormally  thick,  /x3  >>  (*2-  The  dissipation  function  is  given  by 

D  =  P(u,u)  =  \  f\x  K)2  +  fc2(u2)2  +  Hu'tfdz,  =  1  =  (4.3) 

where  ttj',u2  can  be  thought  of  as  the  rates  of  change  of  appropriate  curvatures,  while 
U3  is  the  rate  of  change  of  the  contraction  coefficient  U3.  The  k,'s  are  positive  constants 
reflecting  the  rates  of  energy  dissipation  due  to  deformation  of  material  in  the  beam. 
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Remark  4.1.  It  might  seem  at  a  first  glance  that  the  first  two  terms  in  the  integral 
(4.2)  should  be  replaced  by  a  quadratic  form  (in  u'{,  u") 

aK)2  +  26uy«;  +  e(«5)2; 

however,  by  properly  turning  the  body  coordinate  system  around  the  2-axis,  we  can  di¬ 
agonalize  this  form.  In  Figure  4.2,  the  r-axis  is  chosen  along  the  direction  in  which  the 
beam  bends  most  easily  (pj  <  /i2);  consequently,  the  beam  offers  the  stiffest  resistance  to 
bending  within  the  y2-plane. 


Figure  4.2:  A  beam  with  elastic  coefficients  pj  <  p2. 


Remark  4.2.  We  must  point  out  that  the  beam  was  assumed  to  be  of  uniform 
crossection.  Expression  (4.2)  would  have  to  be  modified  to  include  the  beams  with  variable 
crossections  like  the  ones  show  in  Figure  4.3. 


Figure  4.5:  Beams  with  variable  crossection. 

The  modification  is  in  fact  quite  simple:  one  would  only  have  to  replace  the  first  terms  in 
(4.2)  by  a  quadratic  form  with  r-dependent  coefficients.  In  the  case  of  a  helical  beam  we 


would  take  the  matrix 


(l  c)  =AeM')  =  R-'(‘ S  ®)  * 

where  R(z~J  is  a  rotation  matrix  (  CO?  a  sm  Q  ]  by  2-dependent  angle  a  =  const  •  z 
'  ’  \  —  sin  a  cos  a }  r  ° 

Remark  4.3.  To  incorporate  torsional  deformations  of  the  beam,  we  introduce  the  : 
torsion  angle  a(z\  t),  which  is  the  angle  formed  between  the  i-axis  and  the  projection  onto 
the  xy-plane  of  the  segment  rigidly  connected  to  the  beam  so  that  in  the  unperturbed 
position  of  the  beam  it  is  attached  at  (0,0,  z)T  and  is  parallel  to  the  x-axis.  Thus  the 
variable  o(z:  t)  describes  a  normal  bundle  of  the  beam.  Potential  energy  of  the  beam  is 
given  by 

V  =  \f\<  R~{nRv",v"  >  +p3(u'3)2  +  M4 (a')2)<fc, 
where  p  =  diag(n\,Hl),  M4  is  the  torsional  elasticity  coefficient,  and  t;  =  (u^uz). 

The  point  of  these  remarks  is  that  the  analysis  of  the  rotational  dynamics  of  fairly 
complex  structures  is  possible.  To  keep  our  discussion  maximally  lucid,  however,  we  eschew 
the  complexity  involved  in  more  detailed  structural  mechanical  models,  and  we  refer  to 
the  work  of  Antman  and  Nachman  [2],  for  a  further  analysis  of  the  structural  mechanics 
of  beams. 

The  following  theorem  is  a  straightforward  consequence  of  the  results  presented  in  the 
previous  two  sections. 

Theorem  4.1.  Given  the  system  depicted  in  Figure  4.1  and  described  above  with  kinetic 
energy  (4-1),  potential  energy  (4-2)  and  dissipation  function  (4-2),  the  equations  of  motion 
are  given  by 


Da  +  (m  {,c  +  / 
Jo 

c 

u)  x  Y~xy  =  0 

(4.4) 

D2u  -1-  pdu  +  kdu  +  Y  ly  =  0 

(4.5) 

mb(y  +  Yc)  + 

rt 

/  YD2u  dz  =  0 

'0 

(4.6) 
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where  the  quantities  in  these  equations  are  given  as  follows.  D{  )  —  )  +  w  x  (  )'  a(0  — 

I*>(t)  4-  /o  u  x  ( Du)dz ,  /  »  Me  inertia  tensor  in  the  body  frame  defined  above,  m*,  is  the 
mass  of  the  rigid  body  component,  Y  =  YS~l( w)  (with  S(  •)  is  as  defined  as  in  Lemma 
2.1),  and  d  is  the  differential  operator  defined  by  d  =  [d\ , d\ ,  - d\) ,  p  =  diag(pi,p2,P3) 
and  k  =  zfTffy (fcj ,  k2 ,  k2) . 

Equations  (4.4)-(4.6)  are  equivalent  to 

rt 

Iu  +  u  x  Iu  4-  /  u  x  [ult  +  wxu  +  2wxu(  +  wx(wxu)  +  Y~ly\dz  (4.7) 
Jo 

+mfcc  x  Y~xy  =  0 


un  +  ijxi/  +  2jxu|  +  u;x(uxii)  +  pdu  +  kdut  +  Y~‘y  =  0,  (4.8) 

my  +  Y(J  udz  +  m^c)  =  c'  +  e"t ,  (4.9) 

where  mfc  denotes  the  mass  of  the  rigid  body  component,  and  m  denotes  the  total  mass 
of  the  body-beam  system.  Since  the  mass  density  of  the  beam  has  been  normalized  to  be 
1,  we  find  m  =  t  +  m^.  (4.7)  can  be  further  rewritten,  using  (4.8),  as 

rl 

lu:  4-  u)  x  Iu)  —  J  u  x  \pdu  4-  cdu]dz  4-  m^c  x  Y-1y  =  0.  (4.7)' 

One  easily  checks  that  the  total  angular  momentum  f(y  +  Y  u)  x  j^(y  +  Y u)dm  is  conserved, 
either  using  Noether's  theorem  or  by  direct  computation. 

Remark  4.4.  If  the  center  of  mass  of  the  system  is  at  rest  with  respect  to  the  space  frame, 
we  may  assume  that  c'  =  e"  =  0. 

We  indicate  the  physical  meaning  of  the  terms  in  the  equations  (4.7)-  (4.9).  The  sum 
of  the  first  two  terms  in  (4.7)  is  interpreted  as  the  rate  of  change  of  the  rigid  body’s 
angular  momentum.  The  second  term  in  the  brackets  gives  the  inertial  force  on  the  beam 
due  to  the  body’s  angular  acceleration,  the  third  term  is  the  Coriolis  force,  the  fourth 
is  the  centrifugal  force  and  the  last  term  is  the  D'Alembert  force.  Thus  equation  (4.7) 
expresses  the  conervation  of  the  total  angular  momentum  in  space  expressed  in  the  body’s 


■ 


coordinate  system.  Equation  (4.8)  is  just  Newton’s  law  for  the  beam  expressed  in  the  non- 
inertial  body  frame  (the  D’Alembert  principle) — it  accounts  for  various  inertial  forces. 
Equation  (4.9)  expresses  the  conservation  of  the  linear  momentum  of  the  whole  system.  It 
is  important  at  this  point  to  make  the  following  remark. 

t> 

Remark  4.5  There  is  an  apparent  paradox  associated  with  eq.  (4.7)':  one  might  expect 
to  be  able  to  express  the  integral  term  in  terms  of  u  and  its  derivatives  at  z  =  0,  since  the 
body  feels  the  beam  only  through  the  attachment  point.  As  it  turns  out,  this  expectation  is 
not  met.  In  formulating  our  continuum  mechanical  model  of  the  beam,  we  have  neglected 
certain  effects  such  as  torsional  deformations.  The  implicit  rigidity  in  our  model  leads  to 
this  nonlocal  character  in  the  equation. 


3.5.  Asymptotic  Dynamics  of  a  Rotating  Elastic  Structure 

In  this  section  we  begin  an  analysis  of  the  asymptotic  behavior  of  the  body-beam 
structure  described  in  the  preceding  section.  For  finite  dimensional  dissipative  mechanical 
systems,  LaSalle’s  invariance  principle  [10]  can  be  used  to  show  that  states  asymptotically 
approach  a  minimal  invariant  subset  of  the  zero  set  of  the  dissipation  function  discussed  in 
subection  3.3.  In  subection  3.6,  it  is  shown  that  this  type  of  analysis  may  be  extended  to 
certain  infinite  dimensional  systems  with  features  in  common  with  our  body-beam  model 
described  by  equations  (4.4)-(4.6).  While  we  do  not  prove  that  all  solutions  to  (4.4)-(4.6) 
tend  to  the  zero  set  of  the  dissipation  function  D  as  defined  by  equation  (4.3),  we  do  offer 
a  more  or  less  complete  characterization  of  the  set  of  asymptotic  equilibrium  states  (by 
which  we  mean  the  set  of  solutions  to  (4.4)-(4.6)  for  which  the  equality  D  =  0  also  holds). 
We  prove,  in  particular,  that  in  asymptotic  steady  state  the  beam  displacement  function 
u{ ■  ,  •)  does  not  depend  on  i.  Moreover,  it  is  shown  that  asymptotic  equilibrium  angular 
velocities  are  constant  vectors  parallel  to  the  principal  axes  of  the  steady  state  inertia 
tensor.  This  means  that  the  motion  of  the  system  is  a  pure  rotation  with  no  precession. 
This  is  the  content  of  the  following  theorem. 

Theorem  5.1:  Solutions  o}  (4-7)-(4-9)  which  are  asymptotic  equilibria  (i.e.  solutions 
which  also  satisfy  D  =  0  with  D  as  defined  in  (4-8))  have  the  following  properties: 

(i)  there  is  no  dependence  on  the  time  variable  t  in  the  beam  function:  u  («,0  =  «oo(2).' 

(ii)  the  angular  velocity  u>  is  a  constant 

(Hi)  the  equilibrium  angular  momentum  is  a  constant,  =  Jx u’oo»  tvith  the  equilib¬ 
rium  inertia  tensor  of  the  combined  body-beam  system  given  by 

Jx  —  I  +  f  utuE  -  uuTdz  -  m(C£CmE  -  CmC£)  (5.1) 

Jo 

where  E  =  the  identity  matrix  and  C =  ^(m(,c  +  /q  udz )  «s  the  center  of  mass  of  the 

body-beam  system  (expressed  in  the  body  frame). 

(iv)  equilibrium  rotations  are  aligned  with  a  principal  axis  of  the  equilibrium  inertia 

tensor,  and  thus 


Proof:  To  prove  (i),  note  that  for  asymptotic  equilibria, 

J  fcj(u”)2  +  /t2 ( 2  +  k$(u'3)2dz  =  0.  Hence  u”  =  u'j  =  ii'3  =  0,  and  the  result  follows  as  a 
consequence  of  the  boundary  conditions. 

Using  the  time-independence  of  it,  we  show  that  corresponding  values  of  u>  are  constant. 
Since  ut  =T),  equation  (4.8)  may  be  rewritten 


x  ti  +  u)  x  (u;  x  u)  +  fxdu  +  Y  1  y  =  0. 


Differentiating  with  respect  to  2  at  z  =  0  and  using  the  boundary  condition  ^  =  k  =  ^  0 
at  t  =  0,  we  obtain 

u)xk  +  wx(wxk)  +  ndu2(0)  =  0 
Denoting  the  last  term  by  (ci,C2ic3)^i  rewrite  this  as 


U,’j  —  ^2^3  \  1 

f  C2 

U>2  4-  U>i  W3  1  =  1 

1  _Cl 

w2  +  lj2  J  1 

(5.3) 


Multiplying  the  first  two  components  by  u.'\  and  u2  respectively  and  adding  we  obtain 


C2U)]  —  CjU>2  —  +  U!2U!2. 


(5.4) 


From  the  last  component  of  (5.3)  it  follow’s  that  the  right  hand  side  is  zero,  and  then  (5.4) 
together  with  the  last  component  of  (5.3)  implies  and  u>2  are  constant,  if  Cj  and  c2  are 
not  both  zero. 

To  establish  (ii)  in  the  case  that  Cj  =  c2  =  0,  we  shall  compare  the  constraint 

U/2  ui2  —  C3  (5.o) 

with  the  constraints  provided  by  the  conservation  of  kinetic  energy 

u>T Joc^  =  T  (5.6) 


and  the  conservation  of  magnitude  of  angular  momentum 


If  these  three  polynomial  equations  have  finitely  many  solutions,  then  one  of  them  will  be 
the  desired  time  independent  w*,.  It  thus  remains  to  investigate  the  case  that  the  algebraic 
set  specified  by  these  equations  has  a  component  of  positive  dimension.  To  express  (5.5)- 
(5.7)  in  tefms  of  a  principal  axis  coordinate  system  associated  with  J we  choose  an 
orthogonaTlnatrix  U  such  that  m  =  Uu  satisfies 


aitc’j  4  &i2UJ\u;2  4  Q]3U>iCJ3  4  a2u>2  4  <>23^2^3  +  a3w3  =  c3 


(5.5)' 


Ji-^i  +  h^\  +  =  T 


(5.6)' 


+  >22-’2  +  =  M2 


(5.7)' 


m 


■»0 


« 


Any  algebraic  subset  of  (5.6)',  (5.7)'  having  positive  dimension  must  be  symmetric  with 
respect  to  the  origin.  Hence  we  may  assume  o12  =  aJ3  =  q23  =  0.  Moreover,  because  the 
diagonal  quadratic  forms  (5.5)  and  (5.5)'  are  related  by  a  similarity  transformation,  we 
may  assume  without  loss  of  generality,  that  aj  =  q2  =  1  and  a3  =  0.  If 

1  1  0 

det  ji  j2  j3  ^  0 

j'2  i2  V 2 

h  h  h 

then  each  triple  of  values  e3,T,M2  determines  a  set  of  values  u2,u2  and  uniquely. 
Thus,  as  above,  a;,*,  must  be  constant,  establishing  (ii).  If  the  determinant  is  zero,  either 
Ji  =  J2  =  J3i  >n  which  case  must  be  constant,  or  else  there  must  be  constants  Aj,A2 


such  that 


.*.1 


J}\  h 

^1  h  +  ^2  j 2  ~ 


Neither  Aj  nor  A2  can  be  zero.  The  last  component  implies  that  A]  4  A2j3  =  0.  which  allows 
us  to  eliminate  Aj  from  this  system  of  equations.  Both  j}  and  j2  satisfy  the  polynomial 
equation  -X2jj3  4  A 2j2  -1  =  0.  Either  ji  =  j2  or  else  j}  4  j2  =  j3.  The  second  possibility 
is  ruled  out  by  physical  considerations.  (It  must  always  be  the  case  that  j3  <  j\  4  j2.) 


Hence,  j]  =  j2  and  the  equations  of  motion  are 

-  _  h  ~  h 

U\  —  - : - OJ  2^3 

Jl 

■  _  h  ~  h 

U^2  —  — ~j - 

^3  ~  0 

Since  /3  ^  0,  these  equations  are  not  consistant  with  (5.3)  unless  =  const.  This 
concludes  the  proof  of  (ii). 

(iii)  follows  from  a  direct  computation  using  (4.4)  and  (4.6).  It  may  also  be  obtained 
as  a  direct  consequence  of  the  parallel  axis  theorem. 

(iv)  follows  since  must  satisfy  x  J0 =  0.  | 

Corollary  5.1:  The  asymptotic  equilibrium  beam  function  uoc( •)  and  the  asymptotic 
equilibrium  angular  velocity  j.'x  are  related  by  equations  (5.1), (5. 2)  together  with  the  fourth 
order  system  of  ordinary  differential  equations 

pdu  =  -fi lo(u  -  Cm),  (5.8) 

where  cm  ts  the  center  of  mass  of  the  body-beam  system  in  the  body  frame,  u  =  (i/j ,  t/2,  u3 -f 
2),  and  SUX  =  wx,  and  where  the  boundary  conditions  are  as  prescribed  in  Section  4. 

Proof:  In  light  of  Theorem  5.1,  equation  (4.8)  may  be  rewritten  as 

pdu  =  — u  —  Y'_1y. 

From  equation  (4.9)  (in  which  c'  =  c"  =  0)  we  see  that  y(t)  =  -V^C™.  Moreover,  since 
y  satisfies  Y'(/)  =  Y'^Hoc,  we  have  that  Y'(<)  =  "i oc,n°c f -  Then  y  =  ~y(t)UlcCm  and 
Y'_1y  =  -TllcCm,  proving  the  Corollary. 

Remark  5.1:  (5.8)  prescribes  a  nonlinear  boundary  value  problem  since  na  depends 
on  tioc(-)  through  equations  (5.1)  and  (5.2).  Some  idea  of  the  complexity  involved  in 
explicitly  determining  u00(-)  may  be  gleaned  from  our  solution  to  the  model  problem 
described  by  equations  (6.l)-(6.3)  in  the  following  subsection. 


Remark  5.2:  A  detailed  analysis  of  the  way  in  which  the  dissipative  dynamical  system 
(4. 4)- (4. 6)  evolves  toward  steady  state  is  beyond  the  scope  of  this  paper.  Nevertheless,  a 
heuristic  description  which  makes  contact  with  the  classical  theory  of  rigid  bodies  may  be 
given  as  follows.  For  a  rotating  rigid  body  there  are  two  well  known  conserved  quantities: 
the  kinetitrsenergy  E.  and  the  magnitude  of  angular  momentum  \M\.  If  M  represents  the 
angular  momentum  vector  with  respect  to  the  body  frame,  then  the  energy  E  =  ^ M T I~ 1  Xf 
is  a  quadratic  form  in  M  with  constant  coefficients  and  the  conservation  of  energy  confines 


Figure  5.2:  Momentum  spheres  in  the  conservative 
and  dissipative  cases. 


M  to  an  ellipsoid,  resulting  in  figure  5.2a.  If  a  dissipative  elastic  appendage  is  present 
then  A/,  while  still  confined  to  a  sphere,  will  move  to  smaller  and  smaller  energy  ellipsoids, 
as  indicated  in  figure  5.2b.  For  the  purpose  of  this  heuristic  discussion  we  ignore  the 
infinite-dimensionality  introduced  by  the  beam. 

It  should  be  noted  that  the  basins  of  the  two  sinks  S j  and  S2  in  figure  5.2b  are  inter¬ 
laced.  and  it  is  difficult  to  predict  which  sink  of  XI  will  appear  j/A/(0)  is  near  the  maximum 
energy  point  N .  This  somewhat  delicate  phenomenon  is  said  to  have  been  overlooked  in 
the  design  of  the  Explorer  II  mission.  The  satellite  oriented  itself  along  the  proper  axis, 
but  in  the  direction  opposite  to  the  desired  one. 

It  must  also  be  noted  tnat  the  picture  in  figure  5.2  is  valid  for  moderate  |A/|  only,  and 
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it  undergoes  bifurcations  as  |A/|  is  increased.  These  bifurcations  are  due  to  buckling  of 
the  beam  at  higher  angular  velocities.  Such  a  phenomenon  is  analyzed  in  complete  detail 
in  a  model  problem  in  the  next  subsection. 


3.6.  Bifurcations,  Stability  and  Dissipation  in  a  Model  Problem 


j  In  this  section  we  illustrate  two  interesting  phenomena  —  bifurcations  and  the  trend  to 

the  predetermined  angular  velocities  in  the  simplest  possible  setting.  To  bring  the  under¬ 
lying  phenomena  into  focus,  we  look  at  a  simplified  model,  thereby  minimizing  technical  ’< 

i  .  / 

!  details  while  retaining  several  interesting  features.  We  will  see  that  the  example  consid-  ‘  - 

ered  below  is  an  infinite-dimensional  Morse  system.  It  is  similar  to  the  general  body-beam 

I  system  considered  in  the  preceding  sections.  We  also  point  out  an  unexpected  connec- 

t 

j  tion  between  this  problem  and  that  of  studying  bound  states  of  a  nonlinear  Schrodinger 

equation,  c.f.  [14], 

Figure  6.1  shows  our  model  consisting  of  a  disk  with  a  beam  attached  to  its  center  and 
perpendicular  to  the  disk’s  plane. 

\  *  * 

I  i 

I  i 


Figure  6.1:  A  model  problem  with  1  rotational  d.o.f. 

The  disk  is  constrained  to  rotate  around  the  2-axis  without  friction,  so  that  the  angular 
momentum  of  the  system  is  conserved.  The  beam  is  constrained  to  the  2-u-plane.  and  all 
the  deflections  are  parallel  to  the  u-axis.  Equations  of  motion  of  the  beam,  including 
internal  damping,  are 


pUti  +  puzzzz  "1"  huzzzzf  —  u>  u, 

(6.1) 

u>( J  pu2dz  +  I)  —  M, 

(6.2) 

=  Mo,0  =  u«(i,0  =  u«2(m)  =  o, 

(6.3) 

where  p,  k  charactrize  elasticity  and  damping  respectively,  and  equation  (6.2)  expresses 
conservation  of  angular  momentum,  /  being  the  disk’s  moment  of  inertia.  From  now  on 


27 


we  set  p  —  fj.  =  1  to  simplify  notation.  Steady-state  solutions  will  satisfy  the  ordinary 
differential  equation 

-uzzzz  +  u2u  =  0.  (6.4) 

(6.4)-(6.3)  define  a  nonlinear  eigenvalue  problem  with  M  as  a  parameter,  the  nonlinearity 
lying  in  fhe  u-dependence  of  u>  in  (6.2). 

Remark  6.1.  The  stationary  problem  (6.4)-(6.3)  admits  a  revealing  variational  for¬ 
mulation:  its  solutions  are  the  critical  points  of  total  energy 

E  =  +  J  (u|*  +  w2u  +  u2)dz  (6.5) 

when  the  angular  momentum  M  is  fixed,  as  one  easily  checks.  This  has  a  simple  physical 
explanation:  As  the  beam  vibrates  energy  is  dissipated.  Ultimately  the  virbatioins  are 
damped  out,  and  the  total  energy  of  the  system  is  extremized.  Actually,  for  some  excep¬ 
tional  (in  a  sense  which  we  make  precise  in  Theorem  6.1)  initial  conditions  the  limiting 
energy  value  will  be  not  minimal,  but  rather  critical.  (See  Theorem  6.1  and  Figure  6.4.) 
We  note  also  that 

E  =  ~  jQ  (uZzt)2dz  (6.6) 

Remark  6.2  on  Stability.  The  same  variational  formulation  suggests  a  stability 
criterion  for  the  dynamical  system  (6.l)-(6.3).  Namely,  consider  the  space  S  of  all  triples 
(u,u,u.)  with  angular  momentum  M,  i.e.  satisfying  (6.2).  If  the  stationary  solution  i0  = 
(u.O.w)  minimizes  total  energy  (6.5)  on  S,  then  this  stationary  solution  is  stable,  even 
in  the  undamped  case.  This  means  that  for  all  initial  data  (u,  ut,u>)<=0  in  S  sufficiently 
close  to  the  minimal  (u,0.u;)  (in  the  energy  metric)  the  solution  will  stay  close  to  the 
minimizing  solution  io  for  all  time.  This  is  quite  similar  to  the  standard  minimal  energy 
criterion  of  Lagrange;  here,  however,  the  energy  is  minimal  only  subject  to  the  angular 
momentum  being  constrained.  This  minimality  is  guaranteed  if  the  following  condition 
holds:  Exz  +  XMZZ  >  0  on  the  tangent  space  to  M  =  const,  at  the  point  x0,  where  A  is 
the  Lagrange  multiplier:  Mz  =  A Ez. 


This  is  the  idea  behind  the  Energy-Casimir  method,  which  is  the  second  derivative 
test  in  the  presence  of  the  constraints  and  which  is  applicable  in  the  infinite-dimensional 

28 


i.'.V-V. 


settings.  This  method  has  been  used  to  establish  stability  in  a  variety  of  problems  (see, 
e.g.,  the  work  and  references  in  Krishnaprasad  and  Marsden  [8]). 

A  complete  picture  of  the  behavior  of  our  model  system,  in  particular  of  the  global 
phase  portrait  and  its  bifurcations,  is  given  by  Theorems  6.1  and  6.2  below.  A  perhaps 
surprisingTonsequence  of  this  result  is  that  the  system  selects  a  discrete  set  of  angular 
velocities  for  its  stationary  rotations  independent  of  the  angular  momentum  M!  Another 
aspect  of  this  result  is  that  our  model  is  an  infinite-dimensional  Morse  system. 

To  state  our  first  result  we  combine  (6.1)  and  (6.2)  into  a  single  nonlinear  equation 

+  uZZ2Z  +  kuzzzzt  =  w2(u)u 

where  *v{u)  =  M/(||u||2  4-  /),  ||u||2  =  /0'  u2dz.  This  can  be  rewritten  as  a  system 

*  =  s(“)  =  (-S-,-/[+u!(,),)sfW  M 


Theorem  6.1.  Qualitative  behavior  of  the  system.  Assume  that  k  >  0,  i.e.  the 
damping  is  present. 

There  exists  an  infinite  sequence  0  <  u/j  <  u-i  <...—*  oc  of  preferred  angular  velocities 

associated  with  the  problem  (6.7)  with  boundary  conditions  (6.3)  such  that  if  the  angular 

momentum  M  lies  in  the  interval  (M*,A/*+i)  =  (/u;*,  ),  the  problem  (6. 7), (6.3)  has  k 

distinct  nontrivial  stationary  solutions  (u,t/)  =  (ut(z), 0),  1  <  i  <  k,  and  a  trivial  solution 

( u,v )  =  (0,0),  with  angular  velocities  w(u,)  =  — -  =  u>,,  l  <  i  <  k  and  u;(0)  =  ^f. 

/+Jo  u,dz 

These  solutions  are  depicted  in  Figure  6.2. 


The  solution  with  the  smallest  angular  velocity  is  dynamically  stable,  i.e.  any  solution 
of  (6. 7), (6. 8)  with  nearby  initial  conditions  remains  close  to  it  for  all  time,  in  the  L7 -norm 
given  by 

IIK«)II  =  l!(u1v)ll  =  l  (u2  +  v2)dz. 

Jo 

Stationary  solutions  with  higher  angular  velocities  are  unstable,  and  furthermore,  i-th  so - ! 
lution  (in  the  order  of  increasing  angular  velocity)  ( u,i> )  =  (u,,0)  has  ( i  —  1) -dimensional 
unstable  manifold;  this  dimension  coincides  with  the  number  of  zeroes  of  ut.  In  particular, 
if  M  <  M\  =  I u,’j,  only  one  mode  u  =  0  is  present  and  is  stable.  As  M  crosses  the  bifur¬ 
cation  values  A/,,  the  system  undergoes  a  series  of  pitchfork  bifurcations  shown  in  Figure 


u0  =  0 


Figure  6.2:  Solutions  to  (6.l)-(6.8) 

There  have  been  several  recent  studies  of  the  evolution  of  solutions  to  infinite  dimen¬ 
sional  systems  toward  equilibria  (see  Dafermos  [4]  and  references  therein).  Theorem  6.2 
below  shows  that  under  mild  (physically  reasonable)  assumptions  solutions  of  our  systems 
approach  the  equilibrium  states  described  above. 

Theorem  6.2  Fix  M  >  0  and  let  k  >  0,  say  k  =  1.  Any  solution  w(t)  of  (6.7)  with 
u0(z).u0(z)  £  if^,2)(0'l)  tends  in  the  L2-norm  as  t  —*  oc  to  one  of  the  stationary  solutions 
from  Theorem  6.1  (See  Figure  6.4). 


(WWW 


ms 
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Figure  6.3  5  denotes  stable  equilibria,  and  t/j,U2 _ denote  unstable  equilibria 

with  1,2,...  dimensional  unstable  manifolds.  The  numbers  and 
L2-norm  of  equilibria  depend  on  the  parameter  M. 


Figure  6.4.  Global  phase  portrait  of  the  system  (6. 7), (6. 8).  Indices  show 
the  dimension  of  the  unstable  manifolds .  uhieh  are  depicted  by  arrows. 


Proof:  of  Theorem  6.1. 

1.  Existence  of  k  - f  1  solutions  for  M*  <  M  <  j.  Stationary  solutions  satisfy  the 
nonlinear  boundary  value  problem  (6.3)-(6.4).  There  exists  a  sequence  0  <  uq  <  u>2  <  •  •  • 
of  angular  velocities  for  which  this  problem  has  nontrivial  solutions  (the  eigenfunctions  of 
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with  boundary  conditions  (6.3))  which  we  denote  by  c,et,  c,  arbitrary  and  ||c( (2)  11^2  = 
1.  Define  =  Iuy*  Freedom  of  choice  of  e,  is  used  to  satisfy  the  angular  momentum 
constraint  (6.2),  which  becomes 

u>,(c?  +  J)  =  Af. 

r»- 

This  can  be  solved  for  c,  if  and  only  if  M  >  Iux  =  A/,;  setting  u,  =  c^e,  proves  the  existence 
of  k  solutions  if  M  >  A/*.  One  additional  solution  is  obtained  by  setting  u*+1  =  0,  = 


2.  Stability  of  stationary  modes.  We  restrict  ourselves  to  the  undamped  ( k  =  0)  case, 
as  the  dissipation  does  not  affect  stability  of  stationary  modes.  The  governing  equations 


can  be  written  as 


-  «««  +  (||u||2  +  /)  "  " 


The  dimension  of  the  unstable  manifold  of  the  mode  u,(z)  is  given  by  the  number  of 
positive  eigenvalues  of  the  linearization  of  the  operator  A(u)  at  u  =  u,.  This  linearization 


is  given  by 


/  Af  \2  4 

-*'(«)*’  =  ~a  v  +  (irair^r/J  ■’-no 


^3  (u,  r)tx. 


where  d  =  and  (u,i>)  =  fj  u(z)v(z)dz. 
At  u  =  u,  we  have 


B v  =  A'(u,)v  =  —dAv  -f  ujfv  —  4jJxM  (u,,t/)u, 

The  last  term  in  this  expression  is  a  scalar  multiple  of  the  orthogonal  projection  of  v  onto 
u,.  Comparing  B  with  the  operator  C  given  by 

C  v  —  —  dAv  -f  u;,2v, 

we  see  that  the  spectrum  of  B  can  be  obtained  from  that  of  C  by  replacing  the  eigenvalue 
0  of  C  corresponding  to  the  eigenfunction  u,  by  -4ujJAf_1  ||u,||2.  Thus  the  eigenvalues  of 


B  are 


cj2  -u;?,w2  -u w,2  —  u.,2_ , ,  —  4~’f  A/  ,||u,||2,u;2 


*  M,  is  the  total  angular  momentum  when  the  beam  is  undeflected,  since  the  straight 
beam's  (u  =  0)  contribution  to  the  angular  momentum  is  zero. 


of  which  precisely  i  —  1  are  positive.  This  shows  that  the  dimension  of  the  unstable  manifold 
at  the  stationary  solution  u,  is  (t  —  l). 


Proof  of  Theorem  6.2  consists  of  two  steps: 

(1)  Showing  that  for  any  (  >  0  the  solution  w{t)  enters  and  stays  in  an  (-neighborhood 
(in  L 2  norm)  of  the  zero  set  of  the  dissipation  function.  This  set  is  {(ti,r)  :  /J  vXIdx  = 
0}  =  {(u,v)  :  (u(x),0),  u  6  tf(12)},  see  (6.6). 

(2)  Showing  that  if  w(t)  stays  in  an  (-neighborhood  of  the  zero  set  of  {t>  =  0},  then  u> 
is  (-close  to  one  of  the  stationary  points  of  the  flow,  figure  6.5. 


Figure  6.5  Solutions  which  remain  i-close  to  {u  =  0}  must  be 
%-close  to  a  stationary  point. 

Proof  of  Step  i.  Pick  (  >  0,  and  consider  two  neighborhoods  of  {v  £  0):  >Jt  = 
{(u.f)  :  \\v\\L2  <  (}  and  Nt/'2.  Our  solution  u/(t)  enters  klc/2  for  some  t  >  0  with  necessity: 
otherwise  for  all  t  large  enough  we  would  have 

E  =  -  J  u]ztdz  =  ~  jf  vlzdz<-\0J^  v2dz  <  A0(|)2,  (6.9) 

resulting  in  E(t)  — *  —  oc,  which  is  a  contradiction.  Here  \0  is  the  smallest  eigenvalue  of 
d *  =  [d* I dz*)  with  boundary  conditions  (6.3). 

Thus  to  prove  that  u'(t)  stays  in  >1,  forever  after  some  T  >  0  it  remains  only  to  exclude 
the  possibility  of  infinitely  many  trips  between  /  2  and  the  exterior  ol  r°e  figure  6.5. 

We  will  do  so  by  showing  that  each  trip  results  in  a  loss  of  energy  at  least  A E  >0  - 

depending  only  on  (  and  u;(0)  at  each  crossing,  so  that  u;(t)  can  afford  only  finite  number 
of  trips. 
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Let  f  i ,  #2  be  two  consecutive  crossing  times  of  the  boundaries  of  X,  and  kl(/2.  We  have 


|£<<i)  -£(!,)!  =  I ff  [”Uz |>|(2-<i|a„(^)’,  (6.10) 

where  we  have  used  the  fact  that  /0'  v2dz  >  (|)2  and  the  Poincare  type  estimate  as  in 
(6.9).  It  remains  only  to  provide  a  lower  bound  on  the  trip  time  |/2  —  <j|;  it  is  provided  by 
the  upper  bound  on  the  velocity  w(t)  which  is  the  result  of  the  smoothness  assumption. 

Lemma  6.1.  Any  solution  u;(<)  =  (u(z,t),v(z,t))  of  (6.7)  with  u0(z),t;0(2)  £  //(12)(0, 1) 
has  velocity  bounded  in  the  L2  norm:  i.e.  there  exists  C  =  >  0  such  that  for  all 

t  >  0 

licwor  =  ll»ni,  +  II  -  9'u  -  +  <  c1. 

Proof  of  this  lemma  is  given  in  the  Appendix. 

By  the  choice  of  fj,/2  and  using  Lemma  6.1,  we  obtain 

k  <  |M<2)  -  MMII  =  II  f  f'M*))*!!  <  \h  -  hi  sup  ||F||  <  C\t2  -  fil, 

*  Jh  t> o 

implying 

|t2~ti|  >  2^, 

which  together  with  (6.10)  proves  that  w(t )  stays  in  >1(  for  all  t  >  T( f). 

Proof  of  step  2.  We  show  now  that  if  w(t)  stays  in  N{  for  all  t  >  T{i)  then  it  must 
tend  to  an  equilibrium  point  of  the  flow  (6.7)  -more  precisely,  we  will  show  that  for  every 
6  >  0  there  exists  an  <  >  0  such  that  if  ||v||  <  f  for  all  t  >  T(t)  then  for  some  equilibrium 
solution  wt  =  (ue,0)  we  have  llue  “  ulll»  <  ^  for  all  t  >  T(t).  Our  strategy  is  to  show 
that  if  a  solution  is  r (/)  not  close  to  an  equilibrium,  then  its  velocity  v(t)  must  grow  thus 
taking  it  outside  the  neighborhood  leading  to  a  contradiction. 

Assume  the  contrary:  There  is  some  S  >  0  such  that  for  every  c  >  0  there  is  a 
t0  =  <0(f)  -  ^(f)  sucb  that  although 

IHOH  <  f  for  all  t  >  T(e)  (6.11) 
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we  have 


IMO  -  t*ell  >  <5>  (6.12) 

for  every  equilibrium  ue.  The  idea  is  to  show  that  being  far  from  any  equilibrium  causes 
an  increase  in  velocity  thus  causing  ||v||  >  f  (a  contradiction).  The  details  are  as  follows,  s 
(6.12)  implies  for  /(u)  =  -d4u  +  uj2(u)u: 

||/(u(<0))||  >  A||u(f0)  -  ue||  >  A 6,  (6.13j 

for  some  A  >  0  independent  of  6,  if  we  choose  6  small  enough  (which  we  do). 

We  have:  (cf.  (6.7)) 

HO  -  l'('o)l!i2  >  II  /  /(«(«))^||  -  ||  f  d*v(t)ds\\.  (6.14) 

From  the  proof  of  Lemma  6.1  it  is  clear  that  there  exists  C  >  0  such  that  ||d4f||/2  <  Cf 
for  t  large  enough;  without  loss  of  generality  we  assume  that  this  already  holds  for  t  >  ta. 

We  have  from  (6.14): 

II *’(011  >  II  l‘  f{u(*))ds\\  -  f(l  +  C){t  -  t0)  for  all  t  >  t0  (6.15) 

Jto 

Furthermore,  if  t'.t"  are  sufficiently  large,  and  if 

ll«<>”)  -  «(l')||  <  yfi 

holds,  then 

||/(u(('))  -  /(«(l'))ll  <  jl|/(«(f'))ll 

This  can  be  seen  from  the  proof  of  Lemma  6.1  which  shows  that  all  the  harmonics  of  suffi¬ 
ciently  high  order  decay  exponentially,  with  u’(f)  tending  to  a  finite-dimensional  subspace 
of  L 2  x  L2 .  Applying  this  remark  to  (6.15)  with  t'  =  t01t"  =  f,  we  obtain  the  final  lower 
bound  on  t;,  using  (6.13): 

IMOII  >1(1-  «.)ll/(«(f.))ll  -  <(1  +  C)((  -  <„)  >  (6.16) 

35 


as  long  as 


l«(0  -  «(WII  <  \f(- 


Using  this  estimate  on  the  velocity  v(t),  we  will  show  that  u>(<)  leaves  the  box  ||v||  < 
t,  ||u  —  u(t0)||  <  yfi  through  the  “horizontal”  boundary,  i.e.  ||e(<*)||  =  t  for  some  first  exit 

rr- 

time  <*  >  f0,  see  figure  6.5.  This  will  result  in  the  desired  contradiction. 

To  see  that  II ill  =  (  is  reached  first,  we  note  that  otherwise 


£ 


M<*)  -  «(*o) il  =  yfty 


and 


<(<•  ~‘o)  >||  f  =  MO  -  «(<.)«  =  y/i, 

J  to 

i.e.  the  time  required  to  reach  the  vertical  boundary  is  large: 

1 


Using  this  in  (6.16)— which  is  valid  for  t0  <  t  < 


v/T 


r 


+  -4)11  2  4= 

v« 


2^-«(l  +  C) 


we  obtain 

>  <, 


which  is  a  desired  contradiction,  if  we  pick  <  sufficiently  small.  | 

This  shows  that  given  any  b  >  0  there  exists  an  f  >  0  such  that  if  u'  stays  in  Xlf.  then 
u(z,t)  is  less  than  b  (in  L2[0,  lj)  from  an  equilibrium  solution  ue  of  (6.8).  Since  the  first 
step  in  this  proof  showed  that  for  (-neighborhood  of  the  zero  set  {v  =  0}  of  the  dissipation 
function,  this  proves  Theorem  6.2. 


APPENDIX 


Proof  of  Lemma  6.1. 

Expand  u(x,t)  in  the  orthonormal  basis  {cj(x)}  of  the  eigenfunctions  of  dA  with  the 

boundary  ^conditions  (6.3): 

-  00 

«(*i0  =  H  ai(Oe*(x); 

*  =  o 

system  (6.1)  is  equivalent  to  a  series  of  ODE's  for  the  amplitudes: 


d*  +  =  u>2(1k-. 


{A*}  are  the  eigenvalues  of  dA ,  or 

fo  =  b 

{  b  —  —  Xb  —  Xa  +  w2a, 

where  we  have  dropped  the  subscript  k  for  the  sake  of  brevity. 
It  suffices  to  prove  that  there  exists  C  >  0  such  that 


[A.  1) 


lla'all1  +  =  f;  Al(a|(l)  +  4|(f))  <  C  for  all  I  >  0.  (4) 

k~0 

This  would  show  that  ||F(u;)||£2  is  bounded  for  all  t,  since  the  boundedness  of  the  last 
term  u;2(u)u  in  L2  is  obtained  as  follows: 


u.'2(lt)ll 


(y)4A0-!C. 


To  estimate  a*(<),  &*(0  we  use  the  smoothness  of  the  initial  conditions  to  get  an  upper 
bound  on  a* (0), 6^(0):  estimates 


P12“olli2  =  Ha*(0)Af  <  c 
l|d12v0||  =  £&jfc(0)Af  <  c 


imply 

|a*(0)i,  |6*(0)|  <cAt-3. 


We  will  show  that  for  all  /  >  0  the  estimates 


IMOI  <  CV 


(A.  3) 


hold,  for  k  large  enough,  thus  implying  (A. 2): 


H  A*<»*  <  53  A*c2At  4  =  c2  A;2  <  oo 


since  A*  ~  /r4. 


To  show  (A. 3)  we  rewrite  (A.l)  as 


z  =  Az  +  Rz 


where 


=  (  °  1 
V-A  -A 


M°’  2)- 


Introducing  the  Lyapunov  matrix 


we  obtain 


B  =  [X  eAT‘eAsds  satisfying  AT B  +  BA  = 
Jo 


( Bz.z )  =  -(2,2)  +  ((RtB  -r  BR)z.  z)  < 


■{1-\\RtB  +  BR\\){z,z). 


An  explicit  computation  gives 


1  i  i  _L 

1  +  2*  2A 

J-  —  4-  1 

2A  2 A  T  2A5 


and  thus  ||fZrB  +  BJ?||  <  ,  since  u;  <  y  . 


We  have 


5(B2.2)<-(1- 


“const." 


)(2^)  <  -;>(2>2)- 


the  latter  inequality  holding  for  all  A  =  large  enough. 

Now,  the  minimal  eigenvalue  p  of  the  positive  definite  matrix  B  is  estimated  from 
below  as 

p  >  —  (for  A  large  enough). 

We  have  for  all  /  >  0  and  for  k  large  enough 

<  ( Bz{t),z{t ))  <  ( Bz0,z0 )  <  (20.20), 


al(t)  +  bl(l)  <  2Xk(al(0)  +  b2k(0))  < 

<  2A k2A;6  =  2c2A^5. 


This  impbes  (A. 3),  q.e.d. 


4.  Publications 


The  following  is  a  chronological  list  of  publications  supported  under  this  research  effort. 
All  work  is  jointly  authored  by  John  Baillieul  and  Mark  Levi. 

1.  “The-Mechanics  of  Rotating  Structures,”  MTNS-85:  Seventh  Int’l  Symposium  of  the 
Math.  Theory  of  Networks  and  Systems .  (C.I.  Byrnes  &:  Lindquist,  eds.)  North  Holland 
Publishing. 

2.  “The  Rotational  Mechanics  of  a  Simple  Elastic  Structure,”  Fourth  IF  AC  Symposium 
on  the  COntrol  of  Distributed  Parameter  Systems,  UCLA,  June  30-July  2,  1986. 

3.  “Rotational  Elastic  Mechanics,”  (Preprint). 


5.  Personnel 

Professors  Baillieul  and  Levi  were  the  only  professional  personnel  involved  in  this 
research  this  year. 


6.  Interactions 


Professor  Ballieul  has  given  the  following  lectures  and  presentations  on  this  work: 

1.  “Rotational  Elastic  Dynamics,”  SIAM  Summer  Meeting,  Albany,  NY,  July  17,  1985. 

2.  “The  Invariance  Principle  in  Distributed  Systems,”  IEEE  Conf.  on  Decision  and 
Control  Ft.  Lauderdale,  FL,  December  13,  1985. 

3.  “The  Mechanics  of  Rotating  Structures,”  NSF  Differential  Geometry  Conference,  San 
Antonio,  TX,  April  23-25,  1986. 

During  the  week  of  May  26-30,  1986,  Professor  BaillieuI  will  lecture  on  this  work  during 
his  visit  to  the  Poincare  Institute  at  the  University  of  Grenoble.  Also,  Professor  BaillieuI 
will  give  a  talk  entitled  “The  Rotational  Mechanics  of  a  Simple  Elastic  Structure”  at  the 
Fourth  IFAC  Symposium  on  the  Control  of  Distributed  Parameter  Systems,  UCLA,  June 
30-July  2.  1986. 

Professor  Levi  has  or  will  have  given  the  following  lectures  and  presentations  of  this 
work: 

1.  “Dynamics  of  Flexible  Space  Structures,”  Jet  Propulsion  Laboratory,  Pasadena,  CA, 
December  1985. 

2.  “Dissipation  in  Conservative  Systems,” 

Midwest  Dynamical  Systems  Conference,  Evanston,  11.  April  20,  1986. 

3.  “Bifurcations,  Spacial  Reorientation  and  Dissipation  in  Flexible  Space  Structures,” 
University  of  Maryland,  College  Park,  January  1986. 

4.  “Qualitative  Behavior  of  Flexible  Space  Structures.”  SIAM  Conference  on  Linear  Alge¬ 
bra  in  Signals,  Systems  and  Control.  August  12-14,  1986. 
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